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INFLUENCE  OF  THE  POSITION  OF  AN  ELLIPTICAL  HOLE  ON  STRESS 
CONCENTRATION  IN  A  STRETCHED  ORTHOTROPIC  PLATE 

S.  G.  Lekhnitskiy  and  V.  V.  Soldatov 
(Leningrad,  Ural'sk) 

Problems  in  the  concentration  of  stresses  close  to  holes  in 
anisotropic  plates  are  of  major  practical  interest,  especially 
in  connection  with  the  broad  application  of  synthetic  anisotropic 
materials  in  technology.  They  are  also  of  interest  for  the  minirg 
industry  (study  of  stresses  close  to  underground  shafts,  determination 
of  pressure  on  a  shaft  lining,  etc.).  At  present  approximate 
solutions  of  the  plane  problem  for  anisotropic  media  with  cutouts 
having  widely  varied  configurations  are  known  [1],  but  apparently 
as  yet  there  has  been  no  detailed  investigation  of  the  problem  of 
how  stress  concentration  is  influenced  by  the  position  of  the 
axis  of  the  hole  with  respect  to  the  principal  directions  of 
elasticity  and  the  direction  of  action  of  the  forces.  Exceptions 
to  this  statement  are  the  works  of  A.  S.  Dorogobed  [2]  and 
V.  V.  Lipkin  [3],  where  a  study  was  made  of  the  distribution  of 
stresses  in  a  plate  with  a  round  hole  deformed  by  tangential  and 
bending  forces  which  are  oriented  differently  with  respect  to  the 
principal  directions. 

Given  below  are  some  results  of  an  investigation  of  stresses 
in  an  orthotropic  plate  with  an  elliptical  hole  cut  at  different 
angles  to  the  principal  directions  of  elasticity,  with  one-directional 
stretching. 
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1.  Formulas  for  stress  at  the  edge  of  the  hole.  We  shall 
examine  an  elastic  orthotropic  plate  weakened  by  an  elliptical  hole 
and  stretched  in  one  direction  by  forces  acting  in  the  middle 
plane.  We  shall  make  the  following  assumptions:  1)  the  directions 
x  and  y  of  the  major  axes  of  the  hole  in  general  do  not  coincide 
with  the  principal  directions  of  elasticity  £  and  n;  2)  the  forces 
act  at  an  arbitrary  angle  to  the  principal  axes  of  the  hole,  while 
the  edges  of  the  latter  are  not  fastened  and  are  not  loaded; 

3)  the  material  obeys  the  generalized  Hooke's  law  and  experiences 
small  deformations. 


Fig.  1. 


We  shall  introduce  the  following  designations:  a,  b  -lengths 
of  the  aemiaxes  of  the  ellipse;  c  *>  a/b  (a  >  b);  p  —  intensity  of 
forces;  $  -  their  slope  to  the  major  axis;  -  angle  between 
directions  5  and  x  (Fig.  1);  E.^  and  E2  -  Young's  moduli  for 
tenslon-compre33ion  in  directions  £,  n;  and  v2  -  principal 
Poisson  coefficients;  G  -  shear  modulus  characterizing  the  change 
in  the  angle  between  the  principal  directions 


A"  =- 1'. .m  (1  —  ft)  sin 2\f,  M  -  0.5  [1  —  ft*  +  (1  —  m  +  ft*)  cos  2>pJ  sin  2^ 

'•  =•  o.5n  )1  —  k  +  (1  —  k)  cos  2ip],  N  —  0.25  (1  —  m  V-  k *)  sin*  2^  —  k  (1.2) 

1‘  »  0,125  (3  +  m  +  3**)  +  0.5  (1  —  ft*)  cos  2\ji  4-  0.125  (1  -  m  +  ft*)  costy 

The  equations  of  the  generalized  Hooke's  law  in  the  plane  stressed 
state  are  written  as  follows: 
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1)  for  the  major  system  of  coordinates  q 


2)  for  the  system  x,  y 

*m  “  *11®#  +  *i|#»  +  (nTq 
•*  “  *»l®«  +  «W®#  +  «aTq 

T*»  “  *|A  +  *M®>  +  *«(q 


(1.3) 


(1.4) 


The  coefficients  of  deformation  a^  are  expressed  through  the 
principal  elastic  constants  (see  [1],  pp.  45-47).  In  our 
designations  the  formulas  will  have  the  form 

•«*  +  i)«*2irt 

*»»  *■  l—  v»  +  0.25  (1  —  m  +  Jfc*)  sin*  2ifl 

°-“sr[7+(,-m+*,)«ln*^]  (1.5) 

“  55;  I4  —  4*  +  (t  —m+P) co«2^]>io2^ 

“*•  *  21,  (*  *  ~~  (*  “  *  +  **)  co*  2$]  sin  2$ 


The  equation  of  the  contour  of  the  hole  will  be  taken  in 
parametric  form: 

x-aco«*,  (1.6) 

and  we  will  also  introduce  the  designations 

l*  »  c*  sin’  0  -f  cos’  0 

5’-=  A',[o„c4sii»*»~2rtl,rJ  sin* » con »  +  (2fl„  +  a*) c* sin* 0 008*0—  (1.7) 

—  2a, ,  e  sin  0  con*  +  a„  cos4  0) . 

Assuming  that  the  hole  is  3raall  and  is  located  far  from  the 
edge,  we  will  consider  the  plate  to  be  infinite. 

The  general  solution  for  the  problem  of  stress  distribution  In 
an  infinite  anisotropic  plate  with  an  elliptical  hole  can  be  found 
in  the  book  [1],  In  which  a  number  of  particular  cases  are  also 
examined,  Including  stretching  at  arbitrary  angle  $  ([1],  pp.  148 
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and  152,  formulas  (37.8),  (38.11),  and  (38.12)).  The  expression  for 
the  stress  at  the  edge  of  the  hole  with  arbitrary  anisotropy  is 
quite  complicated  and  is  obtained  only  in  complex  form,  but  in  the 
case  of  an  orthotropic  plate  it  is  possible  to  simplify  it 
significantly  and  to  separate  the  real  part.  We  will  present  the 
final  formula  for  the  stress: 

c,  —  ^  [{CiV  tin*  y  +  (A'  -f  cM)  »lii  ^cos(p  +  (L  +  cP)  con*  q>}  c  sin*  ♦  — 

-  [\MN  +  K(L  +  ei>)l  e tin 1  *  +  [M  (K  -h  cM)  +  cL(L  +  eP)  -  LS)  x 

Xtin?co*i  +  (K  +  cM)Pco*i)l22p£  +  {lK(K+cAf)+  (1,8) 

A'  (N  <—  ££,)]  bIu*  f  •+-  [UN  +  K(L  +  eP)]  aiu  f  cos  ?  +  PN  co»*  j  . 


From  this  (with  $  and  i>  equal  to  0  or  (l/2)ir)  we  will  obtain 
an  expression  for  stresses  for  four  basic  cases,  when  the  hole  is 
cut  In  such  a  way  that  the  direction  of  the  major  and  minor 
axes  of  the  ellipse  coincide  with  the  principal  axes  of  elasticity, 
and  tension  is  carried  out  in  the  direction  of  the  major  or  minor 
axis  of  the  ellipse.  In  particular,  during  stretching  in  the 
direction  of  the  minor  axis  (<p  —  ir/2)  we  will  have 


or 


3*  ^  ■in*  I)  •  i-  (k  +  cn)  cos1 0]  (*“0).  (1.9) 

«. -  [- l.c* »in*a  (1  -i- c/i)co»*0]  (<P~y).  (1.10) 


In  this  case,  when  c  >  1  it  is  natural  to  expect  that  the 
greatest  stress  will  be  obtained  at  the  ends  of  the  major  axis; 
it  is  determined  according  to  one  of  the  two  formulas: 

**  p  (i  +  j j  op  **  0) 

or 

-p(i  -fwi)  (y  «, 

The  stress  on  the  ends  of  the  major  axis  is  determined 
as  compressive  [stress]: 

-  —  /)*:  or  (<>#)►. v,« «  —  £  .  (1  13) 


(1.11) 

(1.12) 
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For  an  isotropic  plate  (k  *  1,  n  *  2)  the  above  assertions 
are  indisputable;  however,  in  the  case  of  an  anisotropic  plate  it 
is  possible  (at  least  in  theory)  to  find  other  relationships  between 
the  constants  at  which  the  greatest  compressive  stress  exceeds  the 
magnitude  of  the  greatest  tensile  stress.  We  will  note  that  stress 
concentration  is  significantly  influenced  by  the  principal  shear 
modulus  G,  entering  into  the  expression  for  n,  which  in  an  orthotropic 
plate  is  not  connected  with  other  elastic  constants;  at  small  G  the 
the  concentration  can  be  very  significant. 

With  an  arbitrary  orientation  of  the  hole  ^  and  stretching  at 
an  arbitrary  angle  of  t  the  theoretical  solution  of  the  problem 
about  which  is  the  greatest  stress  and  where  it  is  obtained  is  made 
difficult  by  the  complexity  of  formula  (1.8).  We  will  limit 
ourselves  below  to  investigation  of  this  problem  for  the  particular 
case  of  a  plate  with  given  elastic  constants. 

All  the  formulas  (1. 8)— (1. 13 )  are  also  valid  for  an  orthotropic 
medium  weakened  by  a  cavity  in  the  form  of  an  elliptical  cylinder 
located  in  a  state  of  plane  deformation  under  the  action  of  stretching 
forces  p;  however  in  this  case  the  moduli  and  E2  and  the 
coefficient  must  everywhere  by  replaced  by  the  corresponding 
quantities : 


*» 

1  1  —  V«V„  • 


Vi'  «=  V,  +  V|,V„ 


(1.14) 


Here  and  v2^  are  Poisson  coefficients  which  characterize 
the  contraction  in  the  direction  of  the  axis  of  the  cavity  during 
stretching  along  the  principal  directions  £  and  n;  and  v^2  are 
coefficients  which  characterize  contraction  along  directions  £  and 
n  during  stretching  along  the  axis  of  the  cavity  (see  [1],  p.  18). 

2.  Results  of  calculations.  At  present  all  of  the  elastic 
constants  which  are  necessary  for  the  calculations  have  been 
determined  only  for  a  comparatively  small  number  of  anisotropic 
materials.  These  include  certain  minerals  (quartz,  topaz,  beryl, 
etc.),  natural  wc  /d  (oak,  pine,  Canadian  fir,  etc.),  delta-wood 
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[veneer  and  plywood],  aviation  plywood,  and  also  SVAM  (glass  fiber 
reinforced  anisotropic  material).  The  Young  moduli  and  Poisson 
coefficients  were  found  experimentally  for  the  latter  for  the  case 
of  stretching  along  the  fibers  and  at  an  angle  of  45  ’  [4],  which 
makes  it  possible  to  calculate  the  principal  shear  modulus;  the 
following  is  obtained: 


'-i  =■  “ ;MC  ' |(t  Kt  i, . '  W*.  C-  n.Ri  ■  10*  (fn  kg/ei^) 

v, ..  v, »  o.  k  ■■■  i  .no,  «  ~ 


(2.1) 


The  following  are  the  greatest  tensile  and  compressive  stresses 
in  a  plate  of  SVAM  with  an  elliptical  hole  in  which  c  -  2  and 
in  which  the  directions  of  the  axes  coincide  with  the  principal 
[directions]  under  tension  along  the  minor  axis  (see 
(1.11M1.13)): 


For  such  an  isotropic  plate  k  -  1,  n  -  2,  a__  _  ■  5p,  and 

UlclX 

o’  _  *  p.  Comparison  of  all  these  numbers  shows  that  SVAM  should 
max 

be  related  to  the  category  of  weakly  anisotropic  materials,  and 
stresses  in  a  plate  of  SVAM  will  differ  insignificantly  from  stresses 
in  an  isotropic  plate. 


Therefore  we  took  a  material  with  stronger  anisotropy,  a  material 
which  has  repeatedly  been  used  to  illustrate  calculations  of 
stresses  in  anisotropic  plates  —  birch  aviation  plywood,  in  which 


El  «=  1.2  •  1C*.  E,  =  0.6  -10*.  C.  =  0.07  •  10*  (in  kg/cn? ) 
v,  =  0.071,  v,=:  0.036,  1.411,  ;i  =  4.453 


(2.3) 


For  the  sake  of  definiteness  we  took  a  =  2b  and  c  =  2.  Tables  1 

and  2  give  the  values  of  the  greatest  tensile  stresses  a  and 

max 

the  greatest  compressive  stresses  o^  for  five  cases  of  tension: 

$  »  0,  30,  45,  60,  and  90c,  and  for  eight  cases  of  orientation  of 
the  hole  with  respect  to  the  principal  directions  of  elasticity: 
ij<  =  0,  30,  45,  60,  90,  120,  135,  and  150°  (accurate  to  2  decimal 
places;  calculations  with  higher  accuracy  are  hardly  advisable, 
keeping  in  mind  the  approximate  character  of  the  initial  data  (2.3)). 
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Table  1.  Greatest  tensile  stresses  (for 

max 

plywood). 


•• 

~i 

M 

!» 

D 

D 

D 

a 

tM 

n» 

0 

3.23 

2.7# 

2.87 

2.21 

2.57 

2.21 

2.87 

2.7# 

3.23 

3.10 

4.82 

4.64 

4.66 

2.55 

3.80 

4.00 

3.76 

3.10 

45 

3.85 

5.S5 

5.84 

6.63 

6.19 

4.28 

4.98 

5.18 

3.85 

Bil 

5.58 

5.47 

6.52 

7.55 

4.24 

5.33 

6.14 

5.58 

80 

7.20 

5.77 

5.52 

7.72 

9.91 

7.72 

5.52 

5.77 

7.30 

Table  2.  Greatest  compressive  stresses  o^,qy/p 
(for  plywood)? 


•* 

♦*-« 

M 

1  ^ 

a 

a 

D 

a 

tM 

1M 

0 

0.71 

0.64 

1.18 

i  .38 

1.41 

1.38 

1.18 

0.64 

0.71 

30 

0.98 

0.77 

0.7S 

1.24 

1.24 

1.54 

1.80 

1.83 

0.08 

45 

1.41 

0.76 

0.46 

1.06 

1.18 

1.38 

1.57 

1.82 

1.41 

60 

1.60 

0.89 

0.63 

1.11 

1.13 

1.24 

1.37 

1.53 

1.60 

90 

1.41 

1.29 

1.13 

0.82 

0.71 

0.62 

1.13 

1.29 

1.41 

We  shall  Introduce  the  maximum  tensile  stresses  for  an  Isotropic 
plate  (for  such  a  plate  In  all  cases  o^ax  *  p): 

**-0  30  45  00  90 

*■«  2.80  IM  4.20  5.00 

Figures  2-5  depict  the  graphs  of  distribution  of  stresses 
o0  along  the  contour  of  a  hole  for  ^  cases: 

1)  9-n>«0;  2)  f  -  90*.  *  -  0;  3)  9- 80*,' 45’;  4)9-30*.  9  =  00*. 

The  magnitudes  of  stresses  were  plotted  from  the  contour  of  the  hole 
on  continuations  of  the  rays  drawn  from  the  center  through  the 
given  points;  positive  stresses  are  depicted  by  arrows  directed  to 
the  center. 

The  diagram  of  loading  is  shown  on  each  graph.  For  comparison, 
a  broken  line  is  used  to  shew  the  graph  for  an  isotropic  plate. 

We  will  note  that  at  one  and  the  same  9  the  general  character  of 
the  graphs  are  retained  if  90°  is  added  to  the  angle  >J>»  i.e.. 
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if  the  axes  C  and  n  change  places,  only  the  magnitude  of  the  stresses 
is  changed. 


Figure  6  Bhows  the  change  in  Cmax/P  (curve  1)  and  °^ax/P 
(curve  2)  as  a  function  of  1>  during  stretching  along  the  minor 
axis  of  the  hole;  Figure  7  gives  the  change  in  omax/p  as  a 
function  of  the  slope  4  of  forces  to  the  major  axis  for  i|>  *  0 
(curve  1)  and  ■  90°  (curve  2). 


iimmi 


By  analyzing  the  tables  and  graphs 
(as  well  as  the  results  of  calculating  stresses 
on  different  points  of  the  contour,  which  we 
did  not  do),  it  is  possible  to  draw  a  number 
of  conclusions  relative  to  a  plate  made  of 
aviation  plywood;  they  give  certain  bases 
for  Judgment  about  the  behavior  under  tension 
of  a  plate  made  from  a  different  material, 
in  which  k  >  1  and  n  >  2.  We  will  note  the 
most  important  of  these. 

1.  The  contour  of  the  hole  is  divided 
into  four  segments  which  are  symmetrical 
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with  respect  to  the  center;  these  segments  are  subjected 
alternately  to  the  action  of  tensile  and  compressive  stresses.  Here 
the  greatest  compressive  stresses  are,  in  all  cases,  significantly 
smaller  than  the  greatest  tensile  stresses. 

2.  In  all  cases  the  maximum  tensile  stress  is  found  to  be 
greater  than  that  in  the  same  isotropic  plate,  or  in  other  words, 
the  concentration  coefficient  for  the  plywood  plate  is  greater  than 
that  for  an  isotropic  plate  (with  equal  +). 

• 

3.  The  greatest  of  all  possible  o  ■  9.91  p  is  obtained  in 
the  case  when  the  opening  is  cut  in  such  a  way  that  the  direction 
of  its  minor  axis  coincides  with  a  direction  of  the  maximum 

Young's  modulus  (fibers  of  the  housing)  and  stretching  is  accomplished 
in  the  direction  of  the  minor  axis. 

4.  Among  the  other  cases  examined,  the  following  turned  out 
to  be  "unsuitable,"  in  the  sense  of  maximum  tensile  stresses: 

1)  s«„»».20/> 

2)  <f-  90*.  oWM  -  7.7 ip 

3)  .  \  a  IKf,  f  s  911°,  Cm,,  “  i.~&p  , 

The  plates  of  greatest  stress  are  located  close  to  the  ends  of 
the  major  axis  (angular  distance  from  the  x-axis  does  not 
exceed  8°). 
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Fig.  6.  Fig.  7. 


5.  During  stretching  along  the  minor  axis  the  least  of 
all  possible  maximum  tensile  stresses,  omax  ■  5.52p,  is  obtained 
in  the  case  when  the  axis  of  the  hole  is  directed  at  an  angle  of 
45°  to  the  principal  directions  of  elasticity. 

Received  28  June  I960 
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*B'Tfi*c(rUNCL,  0)  ABSTRACT  OF  REPORT 

The  problem  of  stress  concentration  near  holes  in  anisotropic 
plates  is  of  considerable  practical  Interest,  especially  in 
connection  witn\  synthetic  anistropic  materials,  with  the  stresses 
gear  undergrounfr/workings .  with  the  determination  of  the  load  on 
“pifc  piKjjJl?';' ^ to’.  The  problem  considered  is  the  stress  around  the 

1  boundary  of  an  elliptical  hole  in  an  orthotropic  plate.  The  plate 

possesses  two  perpendicular  principal  elasticity  directions  in  its 
■  plane;  one  of  these  directions  makes  an  angle  with  the  major  axis 

of  the  elliptical  hole,  and  an  angle  with  the  direction  of  the 
applied  tension,,  jit  is  assumed  that  the  material  obeys  the 
generalized  Hooke »A  law,  that  the  strains  are  small,  and  that  the 
plate  can  be  regarded  as  infinite.  The  formula  for  the  stress 
round  the  boundary/ of  the  hole  is  quoted  from  "Anisotropic  Plates 
by  S,  G,  Lekhnitsmy,  Gostekhizdat,  1957  (Ref.  1),  and  is  discussed 
numerically  using  7the  elastic  constants  of  glass  fibre  anisotropic 
material  and  of  birch  aviation  plywood.  The  glass  fibre  material 
is  weakly  anisotropic  and  the  stress  distributions  do  not  differ 
significantly  from  those  in  an  isotropic  material.  For  the  plywood, 
the  values  of  the  maximum  tensile  and  compressive  stresses  at  the 
edge  of  the  hole  are  calculated,  (The  parent  document  contains 
seven  i graphs  ,  two  tables  showing  greatest  tensile  stresses  and 
greatest  compressive  stresses  for  plywood  which  are  included  on 
microfiche.) 
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